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ABSTRACT 

We describe a comprehensive pulsar monitoring campaign for the Large Area Telescope (LAT) on the Fermi Gamma-ray Space 
Telescope (formerly GLAST). The detection and study of pulsars in gamma rays give insights into the populations of neutron stars 
and supernova rates in the Galaxy, into particle acceleration mechanisms in neutron star magnetospheres, and into the "engines" 
driving pulsar wind nebulae. LAT's unprecedented sensitivity between 20 MeV and 300 GeV together with its 2.4 sr field-of-view 
makes detection of many gamma-ray pulsars likely, justifying the monitoring of over two hundred pulsars with large spin-down 
powers. To search for gamma-ray pulsations from most of these pulsars requires a set of phase-connected timing solutions spanning 
a year or more to properly align the sparse photon arrival times. We describe the choice of pulsars and the instruments involved in the 
campaign. Attention is paid to verifications of the LAT pulsar software, using for example giant radio pulses from the Crab and from 
PSR B1937+21 recorded at Nanjay, and using X-ray data on PSR 10218+4232 from XMM-Newton. We demonstrate accuracy of the 
pulsar phase calculations at the microsecond level*. 
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1. Introduction 

Forty years after the discovery of rotating neutron stars much is 
unknown about their emission processes, and in particular the ra- 
dio emission mechanism is still largely not understood (Lorimer 
& Kramer 2004; Lyubarsky 2008). Of the nearly two thousand 
known pulsars, six have been detected in GeV gamma-rays with 
high confidence, using the EGRET detector on the Compton 
Gamma-Ray Observatory (CGRO) (Thompson et al. 1999). 

The Large Area Telescope (LAT) on the Fermi Gamma-ray 
Space Telescope (formerly the Gamma-ray Large Area Space 
Telescope, or GLAST) went into orbit on 2008 June 1 l(Atwood 
et al. 2008). The sensitivity and time resolution of this instru- 
ment will allow it to discover tens or more of new gamma-ray 
pulsars (Smith & Thompson 2008). Notably, it will be able to de- 
termine the sources among the 169 unidentified EGRET sources 
that are pulsars. However, even with a sensitivity more than 30 



times greater than that of EGRET, the LAT's rate of gamma-ray 
photon detection will be small. For example, the Crab pulsar is 
the third brightest known gamma-ray pulsar, but will trigger the 
LAT only once every 500 revolutions of the neutron star (15 sec- 
onds), on average. While the Crab pulsar should be detected by 
the LAT with high confidence in less than a day, it will take years 
to detect pulsars near the sensitivity threshold, with days sepa- 
rating individual photon arrival times. A search for pulsations 
using gamma-ray data alone is quite difficult in these conditions 
(Atwood et al. 2006; Ransom 2007). Accurate knowledge of the 
rotation parameters increases LAT pulsed sensitivity. However, 
many neutron stars slow down irregularly, a phenomenon known 
as "timing noise", making it difficult to extrapolate a pulsar's ro- 
tation frequency v from one epoch to another Consequently, in 
order to obtain the accurate ephemerides necessary for gamma- 
ray detection of pulsations, known pulsars must be observed reg- 
ularly. 
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In anticipation of the Fermi launch, and mindful of the re- 
quirement for accurate, contemporaneous timing parameters in 
order to observe pulsars at gamma-ray energies, we began an ex- 
tensive campaign of pulsar timing observations with the Parkes 
64-meter radio telescope in Australia (Manchester et al. 2001), 
the Lovell 76-meter telescope at the Jodrell Bank observatory 
near Manchester, England (Morris et al. 2002), and the 94-meter 
(equivalent) Nangay radio telescope near Orleans, France. The 
Parkes telescope is the only telescope in the campaign that ob- 
serves sources south of -39°. Theureau et al. (2005) describes 
the 2002 FORT upgrade to the Nan§ay receiver, with the new 
BON pulsar backend described in Cognard & Theureau (2006) 
and Camilo et al. (2007). These observatories carry out observ- 
ing programs in support of the Fermi mission and, between them, 
observe more than 200 pulsars with a large spin-down lumi- 
nosity, E, as described below, on a regular basis. In addition, 
about 10 pulsars with weak radio emission that are particularly 
strong candidates for gamma-ray emission are being observed 
periodically with the Green Bank radio telescope (GBT) and 
the Arecibo radio telescope. Four pulsars with no detectable ra- 
dio emission are being observed with the Rossi X-ray Timing 
Explorer satellite (RXTE). The Urumqi Observatory (Wang et al. 
2001) is using a 25 meter antenna to monitor 38 of the brighter 
radio pulsars. The goal is to build a database of rotation param- 
eters that will allow folding of the gamma-rays as they are accu- 
mulated over the 5 to 10 year lifetime of the LAT. This work is 
similar in spirit to what was done for CGRO (Arzoumanian et al. 
1994; Johnston et al. 1995; D'Amico et al. 1996; Kaspi 1994). 



2. Pulsars and the Large Area Telescope 

The LAT is described by Atwood et al. (2008). In brief, gamma- 
rays convert to electron-positron pairs in tungsten foil inter- 
leaved with layers of silicon microstrip detectors in the tracker, 
yielding direction information. The particle cascade continues 
in the cesium iodide crystals of the calorimeter, providing en- 
ergy information. Scintillators surrounding the tracker aid rejec- 
tion of the charged cosmic ray background. The scintillators are 
segmented to reduce the "backsplash": a self-veto effect that re- 
duced EGRET'S sensitivity to high energy photons. 

The LAT is a 4-by-4 array of detector "modules" covering 
an area of roughly 1.7 meters on a side. It is sensitive to photons 
with energies between 20 MeV and 300 GeV, whereas EGRET's 
sensitivity fell off significantly above 10 GeV. After event re- 
construction and background rejection, the effective area for 
gamma-rays above 1 GeV is > 8000 cm^ at normal incidence, 
as compared to 1200 cm^ for EGRET. The angular resolution is 
also better than EGRET's, such that source localisation for typi- 
cal sources will be of order of 0.1°' . The height-to-width aspect 
ratio of the LAT is 0.4, for a field-of-view of 2.5 sr, or nearly 
20% of the sky at a given time. Combined with the large effec- 
tive area, this makes a sky survey observation strategy possible: 
on a given orbit, the LAT will sweep the sky 35° away from the 
orbital plane, covering 75% of the sky. At the end of the orbit, 
Fermi will rock to 35° on the other side of the orbital plane, 
and continue to scan. Thus, the entire sky is covered with good 
uniformity every three hours, and no time is lost to earth occul- 
tation. Survey mode, large effective area, and good localisation 
together give the LAT an overall steady point-source sensitivity 
30 times better than EGRET's. 



' Details of the instrument response are maintained at http://www- 
glast.slac.stanford.edu/software/IS/glast_lat_performance.htm 



Gamma-ray events recorded with the LAT have timestamps 
that derive from a GPS clock on the Fermi satelhte. Ground 
tests using cosmic ray muons demonstrated that the LAT mea- 
sures event times with precision relative to UTC significantly 
better than a microsecond (Smith et al. 2006). On orbit, satellite 
telemetry indicates comparable accuracy. The contribution to the 
barycentered time resolution from uncertainty in the LAT's po- 
sition is negligible. 

The EGRET pulsars showed a variety of pulse profiles and 
emission spectra and raised as many questions as they answered 
(Thompson 2004). The high-energy emission is thought to arise 
from basic electtomagnetic interactions of highly relativistic par- 
ticles, namely synchrotron emission, curvature emission and in- 
verse Compton emission. In the two main categories of models 
describing high-energy emission by pulsars, charged particles 
are accelerated along the magnetic dipole field lines by parallel 
electric fields. The "polar cap" model (Ruderman & Sutherland 
1975; Sturrock 1971) argues that the acceleration begins above 
the stellar magnetic pole, but can extend to the outer magneto- 
sphere. In the "outer gap" model (Cheng et al. 1986a,b) particles 
are thought to be accelerated to high energies only in the outer 
magnetosphere, in vacuum gaps between a null-charge surface 
and the light cyhnder. 

The models predict different high-energy emission features 
such as spectra and profiles, that LAT observations may elu- 
cidate, through a hierarchy of observables. First, the different 
models have very different predictions of which and how many 
pulsars emit gamma-rays. Along with detections of radio-quiet 
pulsars in gamma-rays using bUnd search techiuques, the LAT 
analysis using this timing program will constrain the ratio of 
radio-loud to radio-quiet pulsars. This ratio is different for the 
two emission models, with outer gap models predicting a much 
lower ratio (Gonthier et al. 2004; Harding et al. 2007). Reliable 
flux upper limits in the absence of gamma-ray pulsations are use- 
ful in this context (Nel et al. 1996) and also require good timing 
solutions. 

The second observable is the emission profile. Its shape, as 
the beam sweeps the Earth, provides a cross-section of the re- 
gions in the pulsar magnetosphere where the emission origi- 
nates. Coupled with radio intensity and polarization profile stud- 
ies, as well as absolute phase, the gamma-ray light curve pro- 
vides information on the emission geometry, which differs sig- 
nificantly from one model to another (Chiang & Romani 1994; 
Gontiiier et al. 2002). The EGRET pulsars typically have two 
peaks, with the first one slightly offset in phase relative to the 
single radio peak. Although the Crab pulsar breaks this trend, 
LAT observations will study the prevalence of this behaviour as 
a function of pulsar age or other parameters. Pulsar detections 
and emission profiles can only be achieved through solid knowl- 
edge of the pulsar's rotation and good absolute time precision. 
The timing precision will allow finely binned profiles over many 
years even for millisecond pulsars. 

The large energy range covered by the LAT will enable 
measurements of pulsar spectral cut-offs. Although EGRET ob- 
served high-energy cut-offs in pulsar spectra around a few GeV, 
it did not have the sensitivity to measure the exact energy or 
shape of the turnovers. For instance, the LAT should provide a 
determination of the Crab pulsar's spectral cut-off energy, known 
only to be less than a few tens of GeV (de Naurois et al. 2002; 
Teshima 2008), where EGRET lost sensitivity due to the back- 
splash effect. The on-axis LAT energy resolution is better than 
15% above 100 MeV and is better than 10% in the range be- 
tween roughly 500 MeV and 50 GeV, and improves somewhat 
off-axis. The LAT should quickly measure the shape of the Vela 
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pulsar spectral cut-off expected to be around 4 GeV, a power- 
ful discriminator between polar cap and outer gap models and 
a potential diagnostic of high-energy emission altitude (Harding 
2007). Finally, a subset of the pulsars detected by the LAT will 
have sufficient photon numbers to allow phase-resolved spec- 
troscopy, offering further insight into emission mechanisms and 
the beam geometry. 

The LAT will monitor all pulsars continuously with a duty- 
cycle of roughly one-sixth, because of its survey mode, unlike 
EGRET or the Astro-rivelatore Gamma a Immagini LEggero 
(AGILE) telescope (PelUzzoni et al. 2004), which went into or- 
bit in April, 2007. A drawback of the survey strategy is that 
having the sample of gamma-photons spread over a longer du- 
ration makes phase-folding more difficult, as long-term timing 
noise may appear in pulsar spin behavior and glitches may occur 
(Ransom 2007). The need for a substantial and sustained pulsar 
timing campaign stems in part from this continuous observation, 
whereas pointing telescopes only require monitoring during ob- 
servations of any given sky region. 

3. The Timing Campaign 

3.1. Possible Gamma-Ray Pulsars 

For a pulsar with a rotation frequency v (.? ') and frequency 
derivative v = ^ (ii^ units of s^^), the spin-down power is 
E = -An^Ivv erg/s where the moment of inertia / is taken 
to be 10^^ g cm^. The open field-line voltage is V ^ 6.3 x 

lO^'^V^ - 3.18 X 10"^ Ve volts. Above some value of V, 
or, equivalently, E, gamma-ray emitting electron-positron cas- 
cades occur, with ganmia-ray luminosity Ly increasing with E 
(Arons 1996). A linear dependence of Ly on V would give 

leading to a gamma-ray production ef- 
ficiency ey - LyIE oc 1/ Analyses based on EGRET pul- 
sar detections and upper limits have constrained gamma-ray lu- 
minosity laws (e.g. McLaughlin & Cordes 2000), an update of 
which yields Ly oc~ -y-^ '^iP-^, Empirically, although based on 
a small handful of gamma-ray pulsars, the minimum spin-down 
threshold seems to be near £ ^ 3 x 10-'"* erg/s (Thompson et al. 
1999). The angular size and viewing geometry of pulsar beams 
is difficult to constrain and introduces a large uncertainty in the 
relation between a minimum E and the expected gamma-ray 
flux. Bright radio pulsars may have gamma-ray beams missing 
the Earth's line-of-sight; conversely at least one bright gamma- 
ray pulsar, Geminga, has no detectable radio flux (Burderi et al. 
1999). Balancing these issues, and keeping the Ust of gamma-ray 
pulsar candidates of reasonable length, we have selected pulsars 
with E > 10^"* erg/s for LAT pulsar timing. From the ATNF 
online catalogue (Manchester et al. 2005) we obtain 230 such 
pulsars. We give lower priority to the timing of the pulsars in 
globular clusters since they can have apparent E values higher 
than the true spin-down power of the neutron star, due to ac- 
celeration in the gravitational potential of the cluster (Notable 
exceptions to this are the millisecond pulsars PSR B1820 - 30A 
and PSR B 1821 - 24.) This leaves us with 224 pulsars which 
we beUeve are imperative to time regularly. 

Table 1, available electronically, gives the pulsar names as 
well as some indicators of whether they may be gamma-ray 
emitters, such as E and associations with other high-energy 
sources^. The distance d is taken from the ATNF database (the 

^ An up-to-date version is at littps://confluence.slac.stanford.edu/ 
display/GLAMCOG/Pulsars-i-being+timed 



variable "DlSTl"). It is generally based on the NE2001 model 
for the Galactic distribution of free electrons (Cordes & Lazio 
2002) but uses other information such as parallax or HI absorp- 
tion measurements if they are available. The uncertainty in the 
derived distances can exceed 50%, depending on the pulsar. The 

table is sorted by decreasing VS/J^, assuming that oc y as 
discussed above. Such a ranking ignores effects of beam ge- 
ometry relative to the Earth line-of-sight, and variations in Ly 
that may stem from, for example, the angle between the neutron 

star's rotation and magnetic axes. Figure 1 shows nor- 
malized to Vela's value versus the rotation period for the large E 
pulsars. 

Table 1 also hsts some pulsar wind nebulae (PWN) asso- 
ciated with young pulsars (Kaspi et al. 2006; Roberts 2004). Of 
the many striking results recently obtained from the HESS atmo- 
spheric Cherenkov imager array is the large number of Galactic 
sources in the TeV sky, many of which have been identified as 
PWN (see for example Aharonian et al. 2006). Table 1 gives TeV 
associations with HESS sources as well as a MILAGRO source 
(Abdo et al. 2007). Some of the unidentified EGRET sources are 
also hkely to be PWN or pulsars. The table includes the angular 
distances to nearby EGRET catalog sources (Hartman et al. 
1999). Many young pulsars are in or near the error boxes for 
these sources, and the LAT will better localize the GeV sources, 
making coincidence tests stronger. The pulsar timing campaign 
will enhance searches for GeV pulsations, to address whether 
the origin is in the neutron star magnetosphere or in the nebula. 
One study aimed at distinguishing between true and fortuitous 
associations between young pulsars and their PWN or EGRET 
counterparts predicted that 1 9 ± 6 of the EGRET-pulsar proposed 
associations wiU be confirmed by the Fermi LAT observations 
(Kramer et al. 2003). 

The table further lists those rare pulsars seen beyond radio 
wavelengths, either in optical ("O" in the table), or in X-rays. 
The larger gamma-ray pulsar sample expected from the LAT will 
improve the current poor knowledge of the correlations between 
different types of high-energy emission. 

Although we base the LAT timing campaign on high E pul- 
sars, we realize that pulsar gamma-ray emission is far from un- 
derstood and therefore intend to study as many different pulsars 
as possible. The LAT's sensitivity and the continual sky-survey 
mode favor unexpected discoveries. The LAT team therefore 
welcomes long-term, phase-connected rotational ephemerides 
from astronomers wishing to collaborate on pulsed gamma-ray 
searches. 

3.2. Timing Radio-Loud Gamma-Ray Candidates 

The radio telescope time needed to monitor a given pulsar de- 
pends on the precision needed by the LAT, its radio flux density 
(e.g. S 1400 in the ATNF catalog) and pulse profile, and the mag- 
nitude of its timing noise. Simple simulations indicate that gaus- 
sian smearing of gamma-ray arrival times barely degrades de- 
tection sensitivity, for smearing widths up to 0.05 periods. Once 
detected, gamma-photon statistics drive the need for higher pre- 
cision: the timing residuals should be smaller than the phase his- 
togram bin width, which in turn should be wide enough to have 
at least several gamma-photons per bin. 

A consequence of these relatively modest timing require- 
ments is that a given radio observation need only last the min- 
imum time for detection. More crucial is the number of ob- 
servations per year, which depends on the timing noise, corre- 
lated with y and v and thus E (Arzoumanian et al. 1994; Cordes 
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Fig. 1. '^/(f, normalized to Vela, versus neutron star rotation period. 
The pulsars with E > 3 x 10^'* erg/s (squares) and with 10^"* erg/s 
< £ < 3 X lO^'* erg/s (triangles) are those being timed regularly for 
the Large Area Telescope. The 9 solid squares correspond to the six 
confirmed EGRET pulsed detections and three pulsars for which there 
were indications of gamma-ray pulsations. Dots: other pulsars. Pulsars 
in globular clusters are not plotted. Note that radio pulsars with high 
may nevertheless have ganmia-ray beams directed away from 
the Earth line-of-sight and escape detection. 



& Helfand 1980). Gamma-ray candidates tend to be the noisi- 
est pulsars. Illustrations of timing noise in young, high E pul- 
sars can be found, e.g., in Hobbs et al. (2006a). Glitches have 
been observed for roughly a quarter of the pulsars being moni- 
tored (Melatos et al. 2008). The bulk of the pulsars in this cam- 
paign are observed monthly, and a smaller number are observed 
weekly or bi-weekly. 

The low radio fluxes of some gamma-ray pulsar candidates 
require long exposures on the biggest radio telescopes. We must 
devote time to these as they could be bright gamma-ray sources. 
Radio-faint, particularly noisy pulsars could dominate the obser- 
vation schedules. 

Radio signals are dispersed by the interstellar medium, with 
a frequency dependent delay causing signals at high radio fre- 
quencies to arrive before those at low radio frequencies. The 
pulsar Dispersion Measure (DM), or the integrated column den- 
sity of free electrons along the line of sight from a pulsar to 
Earth, usually measured in cm"^ pc, allows extrapolation of the 
photon arrival times from radio to infinite frequency, as is re- 
quired for gamma-ray studies. The DM, however, can change 
over timescales of weeks to years (You et al. 2007). If the 
DM is inaccurate, then the reference phase Oq from the radio 
ephemeris (described below) will change, causing an apparent 
drift in the gamma-ray absolute phase and a smearing of the re- 
sulting gamma-ray pulse profiles. Such smearing would com- 
promise the multi-wavelength phase comparisons upon which 
beam geometry studies are based. Therefore the timing cam- 
paign must include occasional monitoring at multiple radio fre- 
quencies. Figure 2 is one illustration of the magnitude of the 
dispersion for difl'erent radio frequencies. 

Another illustration of the potential effect of DM changes on 
a gamma-ray light curve is obtained using the DM values from 
the Jodrell Bank monthly Crab ephemerides. Over the years of 
the Compton GRO mission (1991-1999), the maximum excur- 
sion in the photon time extrapolated from the radio frequency of 
1400 MHz to infinite frequency is 0.3 ms (1 % of a rotation of the 
neutron star). For reference, the total DM correction from radio 
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Fig. 2. Timing residuals for the 3.05 ms pulsar B1821-24, observed 
with the Nanjay radio telescope using a constant dispersion measure. 
Crosses: 1.4 GHz, Squares: 2 GHz. Adding a dispersion measure time 
derivative to the fit aligns the residuals for the two frequencies. 



to gamma ray is ~ 120 ms, which is 4 Crab rotations. For pulsars 
faster than the Crab, the effect could be larger. For most pulsars, 
the effect is minor. Turbulence in the interstellar medium also in- 
duces frequency-dependent scattering and refraction of the pul- 
sar signal, due to path-length differences. Simulations show that 
those effects are in the order of hundreds of nanoseconds for ob- 
servations at 1.4 GHz (Foster & Cordes 1990), and hence negli- 
gible for gamma-ray astronomy. 

The radio pulsar monitoring must be sustained throughout 
the duration of the Fermi mission (i.e. for 5 to 10 years), a strain 
for any observatory, so other contributions are welcome. In par- 
ticular, very frequent monitoring of high E, large S 1400 pulsars 
could allow significant contributions to LAT science by smaller 
radio telescopes. 

3.3. Radio-Quiet and Radio-Faint Pulsars 

The archetypical radio-quiet gamma-ray pulsar is Geminga, PSR 
J0633-I-1746. Biannual XMM satellite measurements have en- 
sured maintenance of a phase-coherent set of rotation parameters 
over the last few years (Jackson & Halpern 2005). The LAT will 
measure accurate light curves for Geminga in a few days and will 
maintain an accurate ephemeris through gamma-ray timing. The 
AGILE gamma-ray telescope has recently detected Geminga (A. 
PeUizzoni, private communication). 

Table 1 includes at least 15 other pulsars outside of globular 
clusters with S 1400 < 0.1 mJy (some of the others without listed 
S 1400 values are also faint), requiring long radio telescope inte- 
gration times, if detectable at all. 10 of these have E > 10^* erg/s, 
making them both especially promising gamma-ray candidates, 
and subject to especially large timing noise. 

Four high E pulsars are being timed with the Rossi X- 
ray Timing Explorer satellite (RXTE: PSR J181 1-1925 in the 
center of the supernova remnant Gl 1.2-0.3; the young pul- 
sar J1846-0258 in the core of a Crab-like pulsar wind neb- 
ula at the center of the bright shell-type SNR Kes 75, possibly 
the youngest known rotation-powered pulsar (Livingstone et al. 
2006); and B0540-69 and J0537-6910 in the Large Magellanic 
Cloud. 

The remaining six high E, low 5 1400 pulsars are 
J0205+6449, Jl 124-5916, J1747-2958, J1833-1034, 
J1930-H1852, and J2021-H3651. Depending on the pulse 



D. A. Smith et al.: Pulsar Timing for the Fermi Gamma-ray Space Telescope 



5 



shape and the intensity of the surrounding radio nebulae, some 
of these are detectable with the 70-meter class telescopes. For 
the others, the Arecibo and Green Bank (GBT) radio telescopes 
are more appropriate. 

LAT will also perform "blind" searches for new radio-quiet 
pulsars (Ziegler et al. 2008), that is, search for pulsations in 
gamma-ray sources which are not known pulsars. Furthermore, 
the LAT may detect gamma-ray sources bearing neutron star sig- 
natures, for which no pulsations are observed, as was the case 
for the EGRET source 3EG J1835H-5918 (Reimer et al. 2001; 
Halpem et al. 2007). The positional uncertainty obtained with 
the LAT should be small enough so that Arecibo, GBT and 
Parkes can perform deep radio pulsation searches. 

3.4. Public Access to Data 

All pubUc data are released through the Fermi 
Science Support Center (FSSC). The details and 
schedule for the LAT data releases can be found at 
http://fermi.gsfc.nasa.gov/ssc/data/policy/. During the first 
year of operations, which will be devoted primarily to an all-sky 
survey, summary information about a variety of variable sources 
will be released. At the end of this period, all LAT photon data 
and associated science analysis tools will be released; thereafter, 
photon data will be released as soon as processed, typically 
within days of detection. 

As LAT gamma-ray results are published, the ephemerides 
used will be posted on the FSSC server in a "D4 FITS" file (de- 
scribed below). An efibrt will be made to publish a large fraction 
of the timing solutions acquired in this campaign around the end 
of the first year. In any case, the first-year timing solutions for 
all 224 pulsars will be made pubUc 6 months after the end of 
year 1 . Users will be asked to cite the timing parameter creators 
in publications, or to work with them directly. The intent is to 
update a large number of high E pulsar rotation ephemerides in 
the years following, but the continuation of the timing campaign 
will depend on the results of Cycle 1 . 

It is hoped that the pulsar timing data will also be used to 
analyse data from instruments other than the LAT. Atmospheric 
Cherenkov telescopes and neutrino detectors are just two exam- 
ples of experiments that could benefit from the timing campaign. 
Researchers wanting pulsar timing data may contact the authors. 



4. Fermi LAT analysis software 

4.1. The "Science Tools" and f/ie eptiemerides database 

The Fermi LAT "Science Tools" provide a framework for ana- 
lyzing gamma-ray data recorded by the Large Area Telescope: 
data selection, exposure calculation, source detection and iden- 
tification, likelihood analysis of emission spectra, etc The 
"Science Tools" are developed and maintained by the FSSC and 
instrument teams. This software is based on the standard /too/s 
developed at HEASARC ^, designed for data sets using the FITS 
format. In this section we describe gtbary and gtpphase, which 
are pulsar timing analysis tools. 

The pulsar section of the "Science Tools" allows basic timing 
analyses within the FSSC framework, but is not intended to re- 
place specialized packages such as TEMPO (Taylor & Weisberg 
1989) or TEMP02 (Hobbs et al. 2006b). The pulsar science tools 



include only a subset of the functions provided by the those 
packages. 

We have tested gtbary and gtpphase with giant radio pulses 
from the Crab pulsar (B0531+21) and from the millisecond 
pulsar B 1937 +21, recorded at the Nangay radio telescope, X- 
ray photons from the binary millisecond pulsar J0218-H4232 
recorded by XMM, as well as with simulated radio observations 
of the binary millisecond pulsar J0437-4715 from Parkes. 

Furthermore, we made extensive use of pulsar timing solu- 
tions (ephemerides) obtained from radio or X-ray pulsar obser- 
vations. We converted ephemerides to fit the LAT format, a FITS 
file called "D4", which contains a subset of the many parame- 
ters that pulsar astronomers provide. The web interface to the 
ephemerides database generates the D4 FITS file needed by the 
Science Tools. 

When doing a long-term follow-up of a pulsar, one might 
have to use overlapping ephemerides, or choose between 
ephemerides valid on the same epoch. Those ephemerides could 
come from different observatories possibly using different anal- 
ysis methods. As an example, the definition of the arbitrary time 
To when the pulsar rotational phase equals zero, i.e. O(7'o) = 
can differ between observatories. To ensure phase continuity 
when using overlapping pulsar timing solutions, it is important 
to have the template profiles used to build the ephemerides. The 
web-based tool will keep track of these template profiles. In the 
following we describe the pulsar timing analysis using the LAT 
software, and the different tests used to vaUdate this process. 

4.2. Building \ight curves witli LAT software 

Topocentric photon arrival times recorded at the observatory at 
finite frequency have to be transfered to solar system barycenter 
(SSB) times at infinite frequency, mainly by correcting times for 
the motion of the earth and the observatory in the solar system 
frame. Then one folds the barycenter times, using the truncated 
Taylor series expansion for 0(f): 



(5(0 = a>o 



!=0 



fiX(t-Toy- 
(i + iy. 



(1) 



^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/SAE_ovcrvicw.html 
http://heasarc.gsfc.nasa.gov/lheasoft/ftools/ftools_menu.html 



where To is the reference epoch of the pulsar ephemeris, 
fi is the frequency derivative of order i, and Oq is the absolute 
phase, an arbitrary pulsar phase att = To. 

We have tested both barycenter and phase-folding tools. The 
procedure is: 

- Conversion of the arrival times from the observatory-specific 
format to the LAT time format: Mission Elapsed Time 
(MET) TT, which is the number of seconds since 2001 
January 1 at 00:00 (UTC). 

- Calculation of the orbital or ground-based observatory posi- 
tion, and conversion to the LAT spacecraft position format. 

- Transfer of the topocentric times to the barycentric frame, 
using "gtbary". 

- Calculation of the pulsar phase for each arrival time, using 
"gtpphase". 



4.3. Simulated observations of an artificial pulsar 

To test the barycenter software alone, we have simulated arrival 
times at Nangay observatory and compared LAT barycenter soft- 
ware with TEMPO and TEMP 02. 
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Some time and coordinate definition difl'erences exist be- 
tween these different codes. Most TEMPO pulsar timing so- 
lutions have been published using the JPL DE200 planetary 
ephemerides (Standish 1990). TEMPO forms barycentric times 
in "Barycentric Dynamic Time" (TDB). TEMP02 uses the JPL 
DE405 (Standish 1998) solar system ephemerides and com- 
putes barycentric time in "Barycentric Coordinate Time" (TCB) 
units, taking into account the time dilation results from Irwin & 
Fukushima (1999). The LAT bary center tool gtbary handles both 
the DE200 ephemerides and the recommended DE405 model, 
also forming TDB times. The relation between TDB and TCB 
times is given by: 



TDB ^ TCB -LbXAT 



(2) 



where Lb = 1.550519767 x 10"^ ± 2 x 10"" and AT = (date 
- 1977 January 1, 00:00) TAI X86400 s. TAI times refer to 
"International Atomic Time". (TEMP02 has a TEMPO emula- 
tion mode, setting the barycentric time to TDB.) More details 
on time-coordinate definitions can be found in Andersen (1999), 
Rickman (2001) and McCarthy & Petit (2004). 

In the simulation, 10000 arrival times are recorded on the 
ground, beginning on MJD 54100 (arbitrary), with a constant 
step size (no assumption of periodic emission is made), over 1 
year. Nan9ay times are expressed in Modified Julian Days (MJD) 
UTC, at finite frequency. They first have to be moved to the 
LAT time format, at infinite frequency. The dispersion delay in 
the propagation of a signal at a frequency at the solar system 
barycenter fssB through the interstellar medium is the follow- 
ing: 



A?: 



DM 

^fsSB 



(3) 



where K = 2.410 X lO""' MHz"^ cm~^ pc s ' is the dispersion 
constant (see e.g. Manchester & Taylor 1977) and DM is the 
dispersion measure. Note that the frequency at the barycenter 
fssB is difi'erent from the frequency at the observatory, due to 
the Doppler shift resulting from the motion of the observatory 
with respect to the pulsar (Edwards et al. 2006). Higher order 
relativistic corrections are neglected here. The simulated values 
for the pulsar position at J2000 epoch and dispersion measure 
are (a, 6) = (20.75°, 45°), and DM = cm-* pc. 

The position of the radio telescope with respect to the solar 
system barycenter for each time of arrival was calculated us- 
ing the DE200 model in the TEMPO-gtbary comparison, and 
using the DE405 model with TEMP02 in TDB mode for the 
TEMP02-gtbary comparison. The topocentric times are then 
transfered to the SSB. The resulting difl'erences as a function 
of time are shown in Figure 3. In both cases, time differences 
are below 0.7 jus, better than the instrumental precision. We con- 
clude that there is agreement between the LAT barycenter code 
and the other standard tools. 




54100 54150 54200 54250 54300 54350 54400 54450 

Date (MJD) 

Fig. 3. TEMPO-gtbary comparison. TEMPOl-gtbary looks very simi- 
lar. 



detection and times and observatory positions for each date are 
converted to the LAT format. Data were folded using contem- 
poraneous Jodrell Bank ephemerides (Lyne et al. 1993), with 
accuracy better than 160 yus. Event times were then converted 
to the barycenter using gtbary and phase-folded with gtpphase. 
The mean GRP arrival time is 32 yus before that predicted by 
the Jodrell Bank ephemerides, well within ephemeris accuracies. 
The null phase shift of the GRPs relative to the predicted phase is 
consistent with the results of Shearer et al. (2003). This demon- 
strates our ability to phase gamma-ray data over many months, 
even in the presence of significant timing noise, validating the 
codes to a few tens of //s. 

Giant pulses from B1937-H21 were originally discovered 
and studied in detail by Cognard et al. (1996). A study by 
Kinkhabwala & Thorsett (2000) revealed that they occur in win- 
dows shorter than 10 /is, 55 to 70 yUS after the main radio pulse 
and interpulse, allowing us to probe shorter timescales. Dates for 
25 1 giant radio pulses with signal to noise exceeding 30 stan- 
dard deviations were recorded with the Nangay radio telescope 
over three weeks. The timing solution was derived from Nangay 
data. We phase-folded event times having corrected the pulsar 
position for proper motion. Figure 4 shows the resulting phases, 
along with a pulse profile at 1.4 GHz. The mean delays between 
the main and secondary giant pulse components and their regu- 
lar emission counterparts are 60.1 and 67.3 /is respectively, with 
rms deviations of 1 .9 and 2.4 //s, consistent with Kinkhabwala & 
Thorsett's results. The narrow pulse widths demonstrate our pre- 
cision to a few /IS, albeit for a more stable system over a shorter 
duration. 



4.4. Giant radio pulses from the Crab pulsar and B1937+21, 
recorded at Nangay 

Giant pulses are known only from a handful of young and mil- 
lisecond pulsars, and occupy very small windows of pulsar phase 
(see e.g. Johnston & Romani 2004; Knight et al. 2006). Times 
of arrival for 3498 main component Crab GRPs with signal to 
noise ratio exceeding 20 standard deviations were recorded over 
eight months with the Nangay radio telescope, as described in 
Oosterbroek et al. (2008). The radio data were de-dispersed after 



4.5. X-ray data from PSR J021 8+4232, observed by 
XMM-Newton 

Orbital movement has to be taken into account for pulsars in 
binary systems. The 2.3 millisecond pulsar J0218-I-4232 is in a 
binary system with a low mass white dwarf (Kuiper et al. 2002). 
It has been extensively studied at gamma-ray energies and is ex- 
pected to be a bright Fermi source (Guillemot et al. 2007). 

XMM-Newton, an X-ray satellite operating between 0. 1 and 
12 keV, made a 36 ks observation of J0218-I-4232 on 2002 
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Fig. 4. Giant radio pulses from B1937+21 recorded at Nanjay (solid 
line, left-hand scale). The ~ 2 fis pulse width reflects the accuracy of 

the phase-folding. Also shown is the template radio profile at 1.4 GHz Fig. 5. X-ray photons from J0218+4232, recorded by the XMM- 



used for this study (dashed line, right-hand scale). 



Newton satellite, between 1.6 and 4 keV. Black soUd Une: mean back- 
ground level. Black dotted lines: peak positions in Webb et al. (2004), 
respectively Oi = 0.24 and O2 = 0.73. 



February 11-12, with the PN camera''. In timing mode, this in- 
strument has a timing resolution of 30 jUS. Only events well 
calibrated in energy were retained. Pulsar data were collected 
using a rectangular region centered on the source. The back- 
ground level (shown in Figure 5) was estimated by selecting data 
from a similar region, in the same dataset, centered about 50" 
away from the pulsar, where no X-ray source could be detected. 
Finally, only events with energy between 1 .6 and 4 keV were se- 
lected, to allow comparison with studies using the same dataset 
done by Webb et al. (2004). 

Event times recorded by XMM-Newton are expressed in 
MET TT since 1998 January 1 at 00:00 (TT), and hence have to 
be converted to the LAT time format. As for the standard XMM- 
Newton analysis, satellite positions as a function of time were 
determined by a combination of Kepler orbital parameters and 
Chebyshev polynomials. Positions were interpolated to fit the 
LAT position format. 

Event times were converted to the barycenter, then cor- 
rected for the pulsar orbital motion and folded, based on radio 
ephemerides given in Kuiper et al. (2002). We tested frequen- 
cies around the nominal value with a test, and found a shift 
in frequency of Av - 2.6 x lO"*" Hz. As noted by Webb et al. 
(2004), who found a similar offset, such a shift is well within 
the resolution of the timing data. The resulting phase histogram 
between 1.6 and 4 keV is shown in Figure 5. The peaks are cen- 
tered on 0.26 and 0.74 respectively, which are within 40 jUS of 
Webb et al's results. 



4.6. Simulated radio data for PSR J0437-4715 observed 
witti ttie Parkes telescope 

The millisecond pulsar J0437-4715 has a pulse period of 5.76 
ms, and is in a binary system with a 5.74 day orbital period. 
"Post-Keplerian" (PK) parameters, such as the rate of periastron 
advance, d>, or the rate of orbital period decay, Pb, can be fit for 
this binary system (e.g. van Straten et al. 2001). 



5 Wc thank N. Webb (CESR ■ 
XMM-Newton data. 



Toulouse) for providing us with the 



(f498 0.4985 0.499 0.4995 0.5 0.5005 0.501 0.5015 0.502 

Phase 

Fig. 6. Simulated times of arrival from the binary millisecond pulsar 
J0437-4715 recorded at Parkes. The absolute phase is defined so that 
the main radio pulse is centered on 0.5. 



A 500 day observation of J0437-4715 at the Parkes observa- 
tory was simulated using TEMP02 with the FAKE plugin, gen- 
erating 37 times of arrival. The simulation used a timing solu- 
tion for J0437-4715 with 200 ns accuracy, derived from real 
Parkes observations from April 1996 to March 2006 (Verbiest 
et al. 2008). TEMP02 generated a gtpphase-compatihle solu- 
tion, since the LAT Science Tools allow fewer orbital parameters 
than TEMP 02. The accuracy of the simplified timing solution 
was 300 ns. 

Event times were transfered to the solar system barycen- 
ter, having corrected for radio dispersion and pulsar proper mo- 
tion. Times of arrival were phase-folded based on the gtpphase- 
compatible version of the J0437-4715 ephemeris, yielding the 
phase histogram in Figure 6. The mean phase calculated with the 
LAT software is delayed from the TEMP02 mean value by 0.32 
fis, resulting in a vafidation of the code below the /is level. 
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5. Conclusions 

We have motivated and described the large timing campaign that 
is underway for the Fermi mission (formerly GLAST). Previous 
campaigns resulted in a wealth of information on young pulsars, 
and we expect this effort will expand gamma-ray pulsar detec- 
tions to middle-aged, older and millisecond pulsars as well. A 
large database of gamma-ray pulsars of many types will allow a 
study of trends and correlations in important properties such as 
gamma-ray flux, spectral index, profile shape and spectral cutofl". 

Acknowledgements. We made extensive use of the ATNF Pulsar Catalogue 
(Manchester et al. 2005), http://www.atnf.csiro.au/research/pulsar/psrcat/. 
We also made use of the Crab ephemerides provided by the Jodrell Bank 
observatory, http://www.jb.man.ac.uk/pulsar/crab.html. 
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0.56 


4.5e+34 


1.1 


23. 


46. 














J0940-5428 


J0940-5428 


0.50 


1.9e+36 


3.0 


0.66 


— 














J1930+1852 


J 1930+ 1852 


0.44 


1.2e+37 


5.0 


0.06 


— 




X 


3EG J1928+1733 


1.46 




G54.1+0.3 


B1259-63 


J 1302-6350 


0.37 


8.2e+35 


2.8 


1.70 


— 




np 






HESS J1302-638 




J0834-4159 


J0834-4159 


0.36 


9.9e+34 


1.7 


0.19 


— 














J1909-3744 


J1909-3744 


0.36 


2.2e+34 


1.1 


— 


3.* 














B0906-49 


J0908-4913 


0.35 


4.9e+35 


2.5 


10. 


28. 












G270.3-1.0 


J1509-5850 


J1509-5850 


0.34 


5.2e+35 


2.6 


0.15 


— 














B1823-13 


J1826-1334 


0.34 


2.8e+36 


3.9 


2.10 


— 




np 


3EGJ1826-1302(*) 


0.55 


HESS J1825-137 


G18.0-0.7 


J1809-1917 


J1809-1917 


0.34 


1.8e+36 


3.5 


2.50 


— 




np 






HESS J1809-193 




J0538+2817 


J0538+2817 


0.32 


4.9e+34 


1.5 


1.90 


8.20 




X 








S147 


J1811-1925 


J1811-1925 


0.32 


6.4e+36 


5.0 


— 


— 




X 






HESS J1809-193 


Gll.2-0.3, SN386 


J1420-6048 


J 1420-6048 


0.31 


l.Oe+37 


5.6 


0.90 


— 




X 


3EGJ1420-6038 


0.17 


HESS J 1420-607 


G313.6+0.3, Kookaburra 


B 1800-21 


J1803-2137 


0.31 


2.2e+36 


3.9 


7.60 


23. 




np 






HESS J1804-216 




J1046+0304 


J1046+0304 


0.30 


1.4e+34 


1.1 


0.30 


— 












BOl 14+58 


JOl 17+5914 


0.30 


2.2e+35 


2.2 


0.30 


7.60 














J2229+6114 


J2229+6114 


0.28 


2.2e+37 


7.2 


0.25 


— 




X 


3EG J2227+6122 


0.54 




G106.6+3.1 


J1718-3825 


J1718-3825 


0.28 


1.2e+36 


3.6 


1.30 


— 






3EGJ1714-3857 


1.18 


HESS J1718-385 




B 1727-33 


J1730-3350 


0.27 


1.2e+36 


3.5 


3.20 


9.20 






3EG J1734-3232 


1.57 






B0740-28 


J0742-2822 


0.27 


1.4e+35 


2.1 


15. 


296. 














J1617-5055 


J1617-5055 


0.27 


1.6e+37 


6.8 


0.5' 


— 




X 






HESS J1616-508 




J1843-1113 


J1843-1113 


0.27 


6.0e+34 


1.7 


0.10 


— 














J2129-5721 


J2129-5721 


0.26 


2.3e+34 


1.4 


1.40 


14. 














J1124-5916 


J1124-5916 


0.26 


1.2e+37 


6.5 


0.08 


— 




X 








G292.0+1.8, MSH 11-54 


J1846-0258 


J1846-0258 


0.25 


8.1e+36 


6.0 


— 


— 




X 






HESS J1846-029 


G29.7-0.3, Kes 75 


J1913+1011 


J1913+1011 


0.24 


2.9e+36 


4.8 


0.50 


— 










HESS J1912+101 




J1911-1114 


J1911-1114 


0.23 


1.2e+34 


1.2 


0.50 


31. 






3EG J1904-1124 


1.96 






JZU4J+/ /4U 


JZU4J+Z /4U 


0.23 


5.6e+34 


1.8 




1 J. 




X 










J0855-4644 


J0855-4644 


0.22 


l.le+36 


3.9 


0.20 
















J0218+4232 


J0218+4232 


0.22 


2.4e+35 


2.7 


0.90 


35. 




X 


3EG J0222+4253 


1.03 






J 1739-3023 


J 1739-3023 


0.20 


3.0e+35 


2.9 


1.00 








3EGJ1744-3011 (*) 


1.10 






J1831-0952 


J 183 1-0952 


0.20 


l.le+36 


4.0 


0.33 
















B1957+20 


J1959+2048 


0.20 


1.6e+35 


2.5 


0.40 


20. 




np 








G59.2-4.7 
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d 
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% of Vela 


(erg/s) 


(kpc) 


(mjy) 


(mJy) 


Galaxy 


X-ray 


nearby 




assoc. 




Jl 105-6107 


Jl 105-6107 


0.20 


2.5e+36 


5.0 


0.75 






nn 


3EGJ1 102-6103 


0.86 




MSH 1 1-62 


B1821-24 


J 1824-2452 


0.19 


2.2e+36 


4.9 


0.18 


40, 


M28 


X 










B1853+01 


J1856+0113 


0.19 


4.3e+35 


3.3 


0.19 


3,40 






3EG J1856+0114 


0.05 




034.7-0.4, W44, 3C 392 


B 1757-24 


J1801-2451 


0.18 


2.6e+36 


5.2 


0.85 


7,80 




tin 


3EG J1800-2338 


1.26 




G5.27-0.9,G5,4-1,2? 


B0611+22 


J0614+2229 


0.18 


6.2e+34 


2.1 


2.20 


29, 






3EG J0617+2238 


0.69 






B1719-37 


J1722-3712 


0.16 


3.3e+34 


1.9 


3.20 


25, 














J1835-1106 


J1835-1106 


0.16 


1.8e+35 


2.8 


2.20 


30, 














B0540+23 


J0543+2329 


0.15 


4.1e+34 


2,1 


9,00 


29, 














J1913+0904 


J1913+0904 


0.14 


1.6e+35 


3.0 


0,07 
















J0857-4424 


J0857-4424 


0.13 


2.6e+34 


1.9 


0.88 


12, 














J0729-1448 


J0729-1448 


0.13 


2.8e+35 


3.5 


0.70 
















B 1737-30 


J1740-3015 


0.12 


8.2e+34 


2.7 


6.40 


25, 






3EG J 1744-30 11 (*) 


0.86 






J 1928+ 1746 


J1928+1746 


0.11 


1.6e+36 


5.8 


0.25 








3EG J1928+1733 


0.29 






J1015-5719 


J1015-5719 


0.11 


8.3e+35 


5.1 


0.90 








3EG J1014-5705 


0.47 






B 1830-08 


J1833-0827 


0.11 


5.8e+35 


4.7 


3.60 












HESS J1834-087 




J1837-0604 


J 1837-0604 


0.11 


2.0e+36 


6.4 


0.70 








3EG J1837-0606 (*) 


0.18 






J1740-5340A 


J1740-5340A 


0.10 


1.4e+35 


3.4 


1," 




NGC6397 












B 1449-64 


J1453-6413 


0.10 


1.9e+34 


2.1 


14, 


230, 














J 1637-4642 


J1637_4642 


0.099 


6.4e+35 


5.1 


0,78 








3EG J1639-4702 


0,55 






J063 1+1036 


J0631+1036 


0.098 


1.7e+35 


3.7 


0,80 


1,50 




X 










J1702-4310 


J 1702-43 10 


0.096 


6.3e+35 


5.1 


0,72 
















J 130 1-6305 


J 1301 -6305 


0.094 


1.7e+36 


6.7 


0,46 












HESS J1301-631 




J1828~1101 


J1828-1101 


0.092 


1.6e+36 


6.6 


2,90 
















B 1620-26 


J1623-2631 


0.091 


2.0e+34 


2.2 


1,60 


15, 


M4 




3EG J1626-2519 


1,35 






B 1634-45 


J1637-4553 


0.087 


7.5e+34 


3.2 


1,10 


15, 






3EG J 1639-4702 


1,17 






J1702-4128 


J1702-4128 


0.082 


3.4e+35 


4.8 


1,10 












HESS J1702-420 




J1705-3950 


J1705-3950 


0.081 


7.4e+34 


3.3 


1,50 
















J1016-5857 


J1016-5857 


0.080 


2.6e+36 


8.0 


0,46 








3EG J1013-5915 (*) 


0,89 




G284,3-l,8 


J090 1-4624 


J090 1-4624 


0.080 


4.0e+34 


2.8 


0,46 
















B2334+61 


J2337+6151 


0.080 


6.2e+34 


3.1 


1,40 


10, 




no 

"F 










B1317-53 


J 1320-5359 


0.076 


1.7e+34 


2.3 




18, 






3EG J 13 16-5244 


1,75 






J 184 1-0345 


J 184 1-0345 


0.073 


2.7e+35 


4.8 


1,40 








3EG J1837-0423 


1,29 






J1830-0131 


J1830-0131 


0.069 


2.3e+34 


2.6 


0,35 
















Jl 119-6127 


Jll 19-6127 


0.068 


2.3e+36 


8.4 


0.80 






"F 








G292,2-0,5 


J1549-4848 


J1549_4848 


0.068 


2.3e+34 


2.7 


0,47 


17, 














J1816-0755 


Jl 816-0755 


0.065 


2.5e+34 


2.8 


0,17 
















J1857+0143 


J1857+0143 


0.064 


4.5e+35 


5.8 


0,74 








3EG J1856+0114 


0,63 


HESS J1858+020 




BlOOl-47 


J1003-4747 


0.064 


3.0e+34 


2.9 




6,00 














B1610-50 


J 1614-5048 


0.064 


1.6e+36 


7.9 


2,40 
















B1828-11 


J 1830- 1059 


0.060 


3.6e+34 


3.2 


1,40 


2,10 














J1648-4611 


J1648-4611 


0.059 


2.1e+35 


5.0 


0,58 








3EG J1655-4554 


1,68 






B 1643-43 


J1646-4346 


0.057 


3.6e+35 


5.8 


0.98 














G341,2+0,9 


J1906+0746 


J1906+0746 


0.057 


2.7e+35 


5.4 


0,55 


0,90 














BOl 36+57 


J0139+5814 


0.055 


2.1e+34 


2.9 


4.60 


28, 














J1738-2955 


J1738-2955 


0.050 


3.7e+34 


3.5 


0,29 








3EG J1736-2908 (*) 


1,07 






B 1730-37 


J1733-3716 


0.050 


1.5e+34 


2.8 


3,40 
















B0021-72F 


J0024-7204F 


0.049 


1.4e+35 


4.9 


0,15^ 




47Tuc 












J 1723-3659 


J 1723-3659 


0.049 


3.8e+34 


3.5 


1,50 
















J1715-3903 


J 17 15-3903 


0.049 


6.9e+34 


4.1 


0,46 








3EG J1714-3857 


0,33 






J1601-5335 


J1601-5335 


0.048 


l.Oe+35 


4.5 


0,22 
















J1514-5925 


J1514-5925 


0.048 


3.5e+34 


3.5 


0,27 












HESS J1514-591 




B1937+21 


J1939+2134 


0.048 


l.le+36 


8.3 


10, 


240, 




X 










B1820-30A 


J1823-3021A 


0.046 


8.3e+35 


7.9 


0.72 


16. 


NGC6624 












J1828-1057 


J1828-1057 


0.046 


5.5e+34 


4.0 


0.23 
















J0905-5127 


J0905-5127 


0.045 


2.4e+34 


3.3 


1,10 


12, 














Jl 112-6103 


Jll 12-6103 


0.045 


4.5e+36 


12,2 


1,40 
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Notes 


B1221-63 


J 1224-6407 


0.044 


1.9e+34 


3.1 


3.90 


48. 












B 1727-47 


J 173 1-4744 


0.042 


l.le+34 


2.8 


12. 


190. 












B 1356-60 


J1359-6038 


0.041 


1.2e+35 


5.2 


7.60 


105. 












J1921+0812 


J1921+0812 


0.040 


2.3e+34 


3.5 


0.66 














B1915+13 


J1917+1353 


0.039 


3.8e+34 


4.0 


1.90 


43. 












B0021-72E 


J0024-7204E 


0.039 


8.8e+34 


4.9 


0.21'^ 




471^10 










J1838-0549 


J1838-0549 


0.039 


l.Oe+35 


5.1 


0.29 






3EG J1837-0606 {*) 


0.50 


HESS J1841-055 




J 1548-5607 


J1548-5607 


0.039 


8.5e+34 


4.9 


1.00 














J1835-0944 


J1835-0944 


0.039 


5.6e+34 


4.4 


0.41 














J2021+3651 


J2021+3651 


0.038 


3.4e+36 


12.4 


0.10 




X 


3EG J202I+3716 (*) 


0.41 


MGRO J2019+37 


G75.2+0.1 


B1535-56 


J1539-5626 


0.037 


1.3e+34 


3.1 


4.60 














J1638-4608 


J 1638-4608 


0.036 


9.4e+34 


5.2 


0.33 






3EG J 1639-4702 


0.91 






J1910-5959D 


J1910-5959D 


0.036 


5.2e+34 


4.5 


0.05 




NGC6752 










J1841-0524 


J 184 1-0524 


0.035 


l.Oe+35 


5.3 


0.20 






3EG J 1837-0606 (*) 


1.39 


HESS J 184 1-055 




J0024-7204O 


J0024-7204O 


0.033 


6.5e+34 


4.9 


0.10^ 




471^10 










B 1754-24 


J1757-2421 


0.033 


4.0e+34 


4.4 


3.90 


20. 




3EG J1800-2338 


0.97 






J1452-5851 


J1452-5851 


0.032 


3.5e+34 


4.3 


0.24 














J1541-5535 


J1541-5535 


0.032 


l.le+35 


5.7 


0.22 














J1643-4505 


J1643_4505 


0.031 


9.4e+34 


5.6 


0.28 














B 1607-52 


J1611-5209 


0.031 


3.4e+34 


4.3 


1.20 














B1338-62 


J 134 1-6220 


0.030 


1.4e+36 


11.1 


1.90 




X 








G308. 8-0.1 


J1907+0918 


J 1907+09 18 


0.029 


3.2e+35 


7.8 


0.29 


0.40 












J1016-5819 


J1016-5819 


0.029 


4.1e+34 


4.7 


0.31 






3EG J1013-5915 (*) 


1.23 






J0537-6910 


J0537-6910 


0.029 


4.9e+38 


49.4 






LMC X 


3EGJ0533-6916 


1.20 




G279.6-31.7, N157B 


B 1930+22 


J1932+2220 


0.028 


7.5e+35 


9.8 


1.20 


7.80 












J1853-0004 


J1853-0004 


0.028 


2.1e+35 


7.2 


0.87 






3EGJ1856+0114 


1.46 






J0024-7204U 


J0024-7204U 


0.028 


4.6e+34 


4.9 


0.06^ 




471^10 










J1019-5749 


J1019_5749 


0.028 


1.8e+35 


6.9 


0.80 






3EG J1014-5705 (*) 


1.64 






B 1556-57 


J1600-5751 


0.027 


l.le+34 


3.5 


1.40 


20. 












J1841+0130 


J1841+0130 


0.027 


1.2e+34 


3.6 


0.06 














J1543-5459 


J1543-5459 


0.026 


3.8e+34 


4.8 


0.62 














J0954-5430 


J0954-5430 


0.026 


1.6e+34 


3.9 


0.36 














J1737-3137 


J1737-3137 


0.026 


6.0e+34 


5.5 


0.80 






3EG J1734-3232 


1.19 






J1650-4921 


J 1650-4921 


0.025 


1.9e+34 


4.1 


0.16 














J1815-1738 


J1815-1738 


0.025 


3.9e+35 


8.8 


0.25 










HESS J1813-178 




B1822-14 


J1825-1446 


0.025 


4.1e+34 


5.1 


2.60 






3EG J1824-1514 (*) 


0.52 


HESS J1826-148 




B182I-19 


J 1824- 1945 


0.024 


3.0e+34 


4.7 


4.90 


71. 












J 1853+0056 


J1853+0056 


0.024 


4.0e+34 


5.1 


0.21 






3EG J1856+0114 


0.68 






J 1115-6052 


J1115-6052 


0.023 


1.6e+34 


4.1 


0.38 














J1903+0925 


J 1903 +0925 


0.023 


3.2e+34 


4.9 


0.20 














J1413-6141 


JI413-6I41 


0.023 


5.6e+35 


10.1 


0.61 






3EG J 14 10-6 147 


0.80 






J 1756-2225 


J 1756-2225 


0.023 


3.1e+34 


5.0 


0.25 






3EG J1800-2338 


1.50 






J1406-6121 


J1406-6121 


0.022 


2.2e+35 


8.2 


0.36 






3EG J1410-6147 


0.90 






J0024-7204T 


J0024-7204T 


0.021 


2.7e+34 


4.9 






47Hic 










B1718-35 


J1721-3532 


0.021 


4.5e+34 


5.6 


11. 














B 1508-57 


J1512-5759 


0.021 


1.3e+35 


7.3 


6.00 














B 1832-06 


J1835-0643 


0.019 


5.6e+34 


6.2 


1.30 






3EG J 1837-0606 


0.78 






B 1838-04 


J 184 1-0425 


0.019 


3.9e+34 


5.7 


2.60 


2.60 




3EGJ1837-0423 (*) 


1.02 






J1538-5551 


J1538-5551 


0.019 


l.le+35 


7.5 


0.25 














J0024-7204Q 


J0024-7204Q 


0.018 


2.0e+34 


4.9 


0.05.'^ 




47Tuc 










J 11 38-6207 


Jl 138-6207 


0.018 


3.0e+35 


9.7 


0.49 














Jl 156-5707 


J 1156-5707 


0.018 


4.3e+34 


6.0 


0.19 














J1734-3333 


J1734-3333 


0.018 


5.6e+34 


6.5 


0.50 






3EGJ1734-3232 


1.03 






J1412-6145 


J1412-6145 


0.018 


1.2e+35 


7.8 


0.47 






3EGJ1410-6147 


0.53 






J1743-3153 


J1743-3153 


0.017 


5.8e+34 


6.6 


0.50 






3EGJ1744-3011 


1.71 






J1632-4757 


J1632-4757 


0.017 


5.0e+34 


6.4 


0.30 






3EGJ1639-4702 


1.92 


HESS J1632-478 





continuing next page. . . 



continued from previous page 



PSR 


PSRJ 


Ve 


E 


d 


S 1400 


1S400 


Cluster, 


Optical 


EGRET 


EgretDist(°) 


TeV 


Notes 






% of Vela 
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nearby 
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Jl 123-6259 


Jl 123-6259 


0.017 


l.Oe+34 


4.3 


0.56 


11. 














J 1909+09 12 


J 1909+09 12 


0.017 


1.3e+35 


8.2 


0.35 
















J 1650-4502 


J 1650-4502 


0.017 


l.le+34 


4.4 


0.35 








3EG J 1655-4554 


1.41 






J1638-4417 


J1638-4417 


0.016 


3.9e+34 


6.2 


0.21 
















B0540-69 


J0540-6919 


0.016 


1.5e+38 


49.4 


0.02 




LMC 


OX 


3EG J0533-6916 


1.80 




G279.7-31.5, N158A 


J1834-0731 


J 1834-0731 


0.015 


1.7e+34 


5.1 


1.00 








3EG J1837-0606 


1.57 






J 1052-5954 


J1052-5954 


0.015 


1.3e+35 


8.5 


0.15 








3EG J1048-5840 


1.70 






B2148+52 


J2150+5247 


0.015 


l.le+34 


4.6 


2.00 


16. 














J1649_4653 


J1649-4653 


0.015 


l.le+34 


4.7 


0.31 








3EG J1655-4554 


1.71 






J1551-5310 


J1551-5310 


0.015 


8.3e+34 


7.8 


0.54 
















J 1843 -0702 


J 1843-0702 


0.015 


1.2e+34 


4.8 


0.17 








3EG J 1837-0606 


1.85 






J 1907+09 19 


J1907+0919 


0.014 


2.2e+34 


5.8 


















J 1349-6 130 


J1349-6130 


0.014 


1.2e+34 


5.0 


0.58 
















J1913+0832 


J1913+0832 


0.014 


7.4e+34 


7.9 


0.60 
















J1838-0453 


J1838-0453 


0.013 


8.3e+34 


8.1 


0.33 








3EG J1837-0423 (*) 


0.59 






J1627-4706 


J 1627-4706 


0.013 


2.5e+34 


6.1 


0.10 
















J1837-0559 


J1837-0559 


0.013 


1.6e+34 


5.4 


0.50 








3EG J1837-0606 (*) 


0.15 


HESS J1841-055 




J1248-6344 


J 1248-6344 


0.013 


8.6e+34 


8.3 


0.12 
















J 1845 -0743 


J1845-0743 


0.013 


1.3e+34 


5.2 


2.70 
















J 1839-0905 


J1839-0905 


0.012 


1.4e+34 


5.4 


0.16 
















J1845-0316 


J 1 845-03 16 


0.012 


3.9e+34 


7.0 


0.35 












HESS J1846-029 




J1801-2154 


J1801-2154 


0.012 


1.2e+34 


5.3 


0.18 








3EG J1800-2338 


1.75 






J1853+0545 


J1853+0545 


0.012 


1.2e+34 


5.3 


1.60 
















J 1907+0731 


J1907+0731 


0.011 


1.5e+34 


5.7 


0.35 
















B 1557-50 


J 1600-5044 


0.011 


2.8e+34 


6.7 


17. 
















J1755-2534 


J1755-2534 


0.011 


3.5e+34 


7.1 


0.17 
















J1452-6036 


J1452-6036 


0.011 


1.5e+34 


5.8 


1.40 
















J1839-0321 


J 1839-0321 


0.011 


3.6e+34 


7.2 


0.27 








3EG J 1837-0423 


1.22 






J1632-4818 


J 1632-48 18 


0.01 1 


4.8e+34 


7.8 


0.39 












HESS J1632-478 




B 184 1-05 


J 1844-0538 


0.011 


2.3e+34 


6.5 


2.20 








3EGJ1837-0606 


1.83 






J1550-5418 


J1550-5418 


0.011 


l.Oe+35 


9.6 


3.30 
















B2000+32 


J2002+3217 


0.010 


1.2e+34 


5.7 


1.20 


5.50 














J1907+0345 


J1907+0345 


0.0094 


2.3e+34 


7.2 


0.17 
















J1842-0905 


J1842-0905 


0.0091 


l.Oe+34 


5.9 


0.81 
















J1946+2611 


J1946+2611 


0.0089 


l.le+34 


6.1 




1.50 














J1726-3530 


J 1726-3530 


0.0084 


3.5e+34 


8.4 


0.30 
















B 1636-47 


J 1640-47 15 


0.0083 


1.2e+34 


6.5 


1.20 








3EG J1639-4702 


0.38 






J 1904+0800 


J 1904+0800 


0.0082 


3.7e+34 


8.6 


0.36 
















J 1820- 1529 


J 1820- 1529 


0.0080 


4.0e+34 


8.9 


0.61 








3EG J1824-1514 


0.87 






B2011+38 


J20 13+3845 


0.0076 


2.9e+34 


8.4 


6.40 


26. 






3EGJ20 16+3657 


1.95 






J 1908+0909 


J1908+0909 


0.0075 


3.6e+34 


9.0 


0.22 
















B 1855+02 


J1857+0212 


0.0074 


2.2e+34 


8.0 


1.60 








3EG J1856+0114 


1.07 


HESS J1858+020 




J1903+0601 


J1903+0601 


0.0073 


1.5e+34 


7.2 


0.26 








3EG JI903+0550 


0.21 






J1515-5720 


J1515-5720 


0.0072 


l.Oe+34 


6.6 


0.20 
















J1853+0011 


J1853+0011 


0.0070 


2.1e+34 


8.1 


0.30 








3EG J1856+0114 


1.21 






J1342+2822B 


J1342+2822B 


0.0068 


5.4e+34 


10.4 


0.01 




M3 












J 1806-2 125 


J1806-2125 


0.0068 


4.3e+34 


9.8 


1.10 
















J 1626-4807 


J 1626-4807 


0.0066 


2.7e+34 


8.9 


0.37 
















J1020-6026 


J 1020-6026 


0.0065 


9.6e+34 


12.3 


0.14 
















J1305-6203 


J1305-6203 


0.0055 


1.6e+34 


8.5 


0.62 








3EG J1308-6112 


1.08 






J1843-0355 


J1843-0355 


0.0055 


1.8e+34 


8.8 


0.80 








3EGJ1837-0423 


1.59 






J1735-3258 


J1735-3258 


0.0054 


2.4e+34 


9.6 


0.46 








3EGJ1734-3232 


0.66 






B2127+11E 


J2129+1210E 


0.0050 


7.0e+34 


12.9 




0.20 


M15 












B 1758-23 


J1801-2304 


0.0050 


6.2e+34 


12.6 


2.20 








3EGJ1800-2338 (*) 


0.65 






J1043-6116 


J1043-6116 


0.0046 


1.7e+34 


9.5 


0.91 
















J1908+0839 


J1908+0839 


0.0045 


1.5e+34 


9.3 


0.49 

















continuing next page. . . 



continued from previous page 



PSR 


PSRJ 




E 


d 


•51400 


5'400 


Cluster, 


Optical 


EGRET EgretDistf) TeV Notes 






% of Vela 


(erg/s) 


(kpc) 


(mjy) 


(mJy) 


Galaxy 


X-ray 


nearby assoc. 


J1808-2024 


J 1808-2024 


0.0041 


5.0e+34 


13.1 












J1812-1910 


J1812-1910 


0.0035 


1.9e+34 


1 1.2 


0.22 










J 1327-6400 


J 1327-6400 


0.0031 


5.6e+34 


15.5 


0.36 










B2127+11F 


J2129+1210F 


0.0026 


1.9e+34 


12.9 




0.10 


M15 






J 1524-5706 


J 1524-5706 


0.0024 


l.Oe+34 


11.4 


0.41 










J1216-6223 


J1216-6223 


0.0013 


1.3e+34 


16.6 


0.15 










J0535-6935 


J0535-6935 


0.00031 


5.6e+34 


49.4 


0.05 




LMC 




3EG J0533-6916 0.59 


B0456-69 


J0455-6951 


0.00014 


1.2e+34 


49.4 




0.60 


LMC 







Table 1. Pulsars being timed for the Fermi Large Area Telescope (all known pulsars with E > 10^'^ ergs/s), ordered by VS/c/^, where E is the spin-down energy loss rate 

and d is the distance, as an indicator of expected gamma flux suffers many large uncertainties (see text). 5 1400 and 5 400 are the radio flux intensities at 1400 MHz 

and 400 MHz, respectively. "Cluster, Galaxy" is the name of the globular cluster or the host galaxy, if the pulsar is in one. "np" means that the source is observed in X-rays, 
but not pulsed. If the pulsar is located less than 2° away from a 3'''' EGRET catalog source (Hartman et al. 1999), the EGRET name and the angular distance are listed. The 
asterisk (*) indicates that the 3EG source has more than one possible counterpart in the table. The TeV associations are taken from Wagner (2008), available at http://www.mpi- 
hd.mpg.de/hfm/HESS/public/HESSxatalog.htm, and from MILAGRO: (Abdo et al. 2007). Nearby pulsar wind nebulae (PWN) are noted in the last column (Kaspi et al. 2006; 
Roberts 2004). The data in the first 8 columns were obtained from the ATNF database except for the radio flux densities with the superscripts a: F. Camilo, private communication, 
b: Jacoby et al. (2003), c: Kaspi et al. (1998), d: Ray et al. (1996), e: D'Amico et al. (2001),/: Camilo et al. (2000). This table is available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-HA/ 



